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Hot subdwarfs are evolved low-mass stars that have survived core helium ignition and are now in (or recently finished 
with) the core helium burning stage. At the hot end of the Horizontal Branch (HB), many of these stars are multiperiodic 
pulsators. These pulsations have revealed details of their global and internal structure, and provide important constraints 
on the origin of hot HB stars. While many features of their structure deduced from seismic fits have confirmed what we 
expected from evolutionary considerations, there have been some surprises as well. 



1 Introduction 

The fundamental picture of low-mass stellar evolution in- 
cludes the still-enigmatic transition from the red giant phase 
to the hehum core burning phase. Stars that begin life with 
masses less than about 2.3 Mq enter this phase through ig- 
nition of helium while the core is degenerate - the core he- 
lium flash (first studied by Mestel 1952 and Schwarzschild 
& Harm 1962). Core helium-burning models that have sur- 
vived the flash populate the Horizontal Branch (HB), as as 
noted first by Hoyle & Schwarzschild (1955). 

Historically, computation of evolution through the core 
helium flash has been a difficult problem. To get around the 
core helium flash problem, HB models are generally con- 
structed as "Zero Age Horizontal Branch" (ZAHB) mod- 
els using a pure helium core surrounded by a hydrogen-rich 
envelope. To account for the distribution of core helium- 
burning stars in the H-R diagram, the stars must lose a 
significant amount of mass at (or before) the onset of the 
flash, as was first demonstrated numerically in the last 1960s 
(Faulkner & Iben 1966, Rood 1970, Iben & Rood 1970). 
Families of ZAHB models share the same core (the mass 
determined by the size of the hydrogen-exhausted core at 
the onset of the helium flash, with the amount of remnant 
envelope is treated as a free parameter The thinner the hy- 
drogen surface layer, the higher the effective temperature. 

This approach to modeling HB stars has remain largely 
unchanged over the past 40 years. This way of making mod- 
els of stars avoids the helium core flash questions, yet pro- 
vides satisfactory (though not complete) models of HB stars 
in globular clusters and other old stellar populations. How- 
ever, we learn little about the conditions of the core helium 
flash or about the mass loss event(s) that reduce the mass of 
the envelope when modeling stars in this traditional way. 

What has changed over the past 10-15 years is the in- 
creased interest in the origin and evolution of the hottest 
horizontal branch (EHB) stars, such as the hot subdwarf B 
(sdB) stars. To account for the bluest EHB stars, the rem- 
nant hydrogen layer must be thinner than about 0.5% of the 
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Star's mass. Stars with envelopes that thin (or thinner) do 
not support hydrogen shell burning - the entire luminosity 
of the star is generated by helium core burning. Figuring 
out how a star can leave itself with such a thin hydrogen 
layer has perplexed astronomers for decades. As far back as 
1976, Mengel et al. (1976) proposed binary mass exchange 
as one possible mechanism; other proposals for their ori- 
gin include variations of mass loss efficiency on the RGB 
(D'Cruz et al. 1996), and even mass stripping by plane- 
tary companions (Soker & Harpaz 2000). Growing obser- 
vational evidence now suggests that binarity plays a role in 
the formation of the subdwarf B (sdB) stars (Maxted et al. 
2001, Heber 2009). Theoretical support for this view comes 
from a series of simulations of binary evolution (Sandquist 
et al. 2000, Han et al. 2002, 2003). 

Over this same time, discovery that some sdB stars un- 
dergo nonradial pulsation has brought the tools of astero- 
seismology to the study of these stars. Two broad classes are 
known - short-period p-mode pulsators, and longer-period 
g-mode pulsators. and the frequencies provide a probe of 
their global properties and internal properties. These anal- 
yses have concentrated on as well as structural parameters 
such as the thickness of the surface hydrogen-rich layer via 
matching the theoretical p-mode frequencies of models to 
the observations. There are several excellent recent reviews 
of this exciting field - see, for example, Charpinet et al. 
(2008a) and 0stensen (2010). For a general review of ob- 
servations HB stars, see the excellent review by Catelan 
(2009), and for a recent account of all things hot subddwarf, 
see Heber (2009). 

In this brief review, I recount the basic theoretical as- 
pects of sdB interior evolution, and focus on how asteroseis- 
mology addresses questions about the origin of sdB stars, 
and more generally, the physical processes that operate in 
the cores of HB stars. Seismic analysis of p-modes has re- 
vealed the envelope structure of many sdB stars. But some 
of these p-mode pulsators, along with a growing number of 
longer period g-mode pulsators, allow us to pursue the goal 
of probing the cores of sdB stars. In that way we can "look" 
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Fig. 1 Anatomy of a representative sdB evolutionary track 
for a model with a mass of O.478M0. Evolutionary stages 
and the time following the Zero Age EHB are indicated. 
At the start of semiconvection, the core helium abundance 
has been reduced to 0.66. When the model reaches the min- 
imum log 5, that abundance has dropped to 0.12. At that 
point the model has a fully convective core out to a mass 
fraction of 0.28 M, and the top of the semiconvective re- 
gion is at a mass fraction of 0.5 M*. 



at the region of the star that was directly involved in the 
core helium flash. The nonradial pulsation spectrum in this 
subgroup of stars can help us test these models and see the 
physics in action. 



2 Models of sdB evolution 

Ordinary HB stars (those that have envelopes thick enough 
to support hydrogen-shell burning) can arrive on the ZAHB 
through single star evolution in a relatively simple way, such 
as through steady mass loss on the RGB, or a single episode 
of modest mass loss at he core helium flash. Arrival on the 
EHB with very thin envelopes is more problematic. Fine- 
tuning of the mass loss mechanism is needed to produce 
such a thin hydrogen-rich envelope. However, single-star 
scenarios that result in EHB stars produce structures that 
are entirely analogous to those with thicker envelopes. 

2.1 Standard model(s) 

For standard sdB models - with hydrogen-rich envelopes 
that are too thin to support a hydrogen burning shell - the 
evolutionary tracks have a distinctive shape. A representa- 
tive example is shown in Fig[T] constructed as described 
in Kawaler & Hostler (2005). Increasing (decreasing) the 




Fig. 2 Internal structure of two sdB models. Top: helium 
abundance (solid line) and normalized energy generation 
rate per gram (dashed line) for the 96 Myr model that is 
nearing core helium exhaustion. Bottom: the same quanti- 
ties during the shell narrowing phase - the energy genera- 
tion rate peaks in the narrowing helium-burning region sur- 
rounding the C/0 core. The dotted line in the bottom panel 
shows the composition profile of the upper model for com- 
parison. 



thickness of the surface hydrogen layer translates the curve 
towards lower (higher) temperatures and gravities. 

Upon settling on the ZAHB, the core has commenced 
burning helium via the 3q; process. The basic time scale 
for core helium burning in these stars is about 10^ years. 
This highly temperature-sensitive reaction results in vigor- 
ous core convection that generates a steep composition gra- 
dient at the edge of the convective core. One consequence of 
the steep composition gradient is that the convection bound- 
ary is unstable to small amounts of outward mixing. Thus, 
by 30 Myr after core helium ignition, the core is likely to 
overshoot, resulting in extra mixing. If it occurs gently, this 
extra mixing resembles diffusive mixing, where the com- 
position gradient is "softened" to ensure neutral stability 
against convection - a process often called semiconvection. 
The precise mechanism by which this extra mixing occurs 
is not fully understood, even if it can be computed conve- 
niently through mixing to convective neutrality. 

Through most of helium core burning, the effective tem- 
perature of the model stays roughly constant, but the star 
grows in radius and therefore moves upwards in the H-R 
diagram (see Fig.[rii. As the core becomes significantly de- 
pleted in helium, the evolutionary track turns around, and 
the model moves more quickly to higher temperature and 
gravity. The luminosity and helium abundance profile at the 



© 0000 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-joumal.org 



Astron. Nachr. / AN (0000) 



3 



turning point are shown in the upper panel of Fig. |2] At this 
point, the fully convective core reaches out to a mass frac- 
tion of 0.3M, (or O.I4M0), while the semiconvective zone 
reaches the midpoint in mass (0.5M, or O.24M0). 

Following helium exhaustion in the core, a thick helium- 
burning shell develops, working its way quickly out from 
the center as it consumes the surrounding, helium-depleted 
material. As the shell encounters the helium-rich material, 
it slows, and the luminosity (and radius) begin to increase 
(the shell-stalling phase in Fig.[TJ and the model makes its 
way, eventually, to the white dwarf cooling track. 

2.2 Evolution from binary progenitors 

As mentioned above, a leading candidate for the production 
mechanism for sdB stars involves common-envelope evo- 
lution, leading to mass transfer and stripping of the bulk of 
the hydrogen-rich envelope. Stars produced via this mech- 
anism can arrive in the sdB region of the H-R diagram with 
total masses and envelope masses that are significantly dif- 
ferent than for single-star production scenarios (Han et al. 
2002, 2003, Hu et al. 2008). 

Hu et al. (2008) explore evolution to and beyond the ZA- 
EHB using higher-mass initial models that have had mass 
stripped while on the RGB (as a consequence of binary evo- 
lution). Their initial mass is sufficiently high that they ignite 
helium non-degenerately, avoiding the helium core flash but 
still evolving to the hot horizontal branch. Their subsequent 
evolution includes an active hydrogen-burning shell, and 
their evolutionary path is quite different than the standard 
model described above (see Hu et al. 2008 for details). 

2.3 Diffusion 

Because of their relatively high surface gravities and ra- 
diative envelopes, sdB stars are subject to diffusive proper- 
ties that can significantly alter their surface abundances and 
subsurface composition profiles. Observational evidence in- 
cludes significant helium depletion in EHB stars by a factor 
of 10 to 100 (Heber et al. 1986) which likely results from 
gravitational settling. Theoretical calculations can explain 
the general trends suspected for diffusion (depletion of he- 
lium and some heavier elements) illustrating the importance 
of diffusion on the evolutionary paths of these stars for both 
standard models and those with binary progenitors (i.e. Hu 
et al. 2010). 

In addition to gravitational settling, another important 
process is radiative levitation of metallic species (Fontaine 
& Chayer 1997, Moehleret al. 2000, Unglaub & Bues 2001), 
which can account for peculiar distribution of these ele- 
ments in the spectra of these stars (Geier et al. 2010). Inclu- 
sion of even very small amounts of stellar wind can also al- 
ter the surface abundances (Fontaine & Chayer 1997, Unglaub 
& Bues 2001, Chayer et al. 2004, Unglaub 2008) and in fact 
may be necessary to understanding the detailed abundance 
patterns. Models that include all of these diffusive processes 
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also predict that iron-peak elements should also accumulate 
below the surface of these stars which turns out to be crucial 
to explaining the driving mechanism for the pulsating sdB 
stars (Charpinet et al. 1997). 

3 Seismology of sdB stars - envelopes 

Our desire to understand the origin of sdB stars requires 
more information than can be provided by traditional spec- 
troscopic observation. Such observations yield precise posi- 
tions in the H-R diagram which can be compared, both indi- 
vidually and as an ensemble, to evolutionary tracks. Time- 
resolved spectroscopy and photometry can also reveal the 
presence (and properties) of binary companions as well. But 
it is hard to construct a detailed make-or-break experiment 
that could be decisive in confirming a production scenario 
that is based on this approach. 

However, as shown by Hu et al. (2008), the internal 
structure of an sdB star can reveal its progenitor by its in- 
fluence on the observed p-mode spectrum. We can, for ex- 
ample, use asteroseismic probes will be able to measure the 
position below the surface of the H/He transition layer for 
comparison to models with different origins. 

In addition, these probes also allow us to test models of 
diffusion, radiative levitation, and winds in sdB stars. Com- 
bined with spectroscopic measurements of elemental abun- 
dances, this powerful suite of observational tools can give 
us an empirical picture of the physics operating within the 
envelopes of sdB stars. 

3.1 Hydrogen layer thickness - adiabatic 
asteroseismology 

The short-period sdB stars are pulsating in p-modes, which 
sample the surface layers much more sensitively than the 
core. Therefore, comparison between observed pulsation fre- 
quencies and models (constrained by spectroscopically de- 
termined values of Tcff and log g) can reveal details about 
the surface layers (notably, Mh) as well as the global prop- 
erties (mass, radius, and luminosity) of the stars. 

The most exhaustive efforts at seismic probing of the 
short period sdB pulsators is the work over the past decade 
by Charpinet, Fontaine, Brassard, and their collaborators. 
They have analyzed seismic data on more than a dozen of 
the short-period sdB stars, and determined estimates of the 
total mass and surface hydrogen layer thicknesses for most 
of those stars. A good recent summary of the many results 
from their analysis can be found in Fontaine et al. (2008); 
see Ostensen (2010) for an updated summary plot showing 
the distribution of resulting core/envelope masses. 

The determined range of masses is wider than would be 
expected for single star evolution, suggesting that the bi- 
nary formation channel is indeed an important contributor 
to the observed sdB population. While further work needs 
to be done to estimate the relative contribution of the binary 
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channel(s) this result demonstrates the potential of astero- 
seismology to address this fundamental question about the 
origin of sdB stars. An additional important conclusion is 
that the envelope thickness determinations are, in all cases, 
as small as the standard model suggests - indicating that, 
as expected, the surface hydrogen layer in these stars is too 
thin to support a hydrogen burning shell. 

3.2 Diffusion and "nonadiabatic" asteroseismology 

Charpinet et al. (1997) were the first to show the importance 
of subsurface iron in driving pulsations of sdB stars (see also 
Fontaine et al. 2006, Charpinet et al. 2008a, and references 
therein). Without radiative levitation, the iron abundance is 
too low for the iron-peak opacity to drive the oscillations. 
Thus the very existence of pulsations in sdB stars is evi- 
dence that radiative levitation of heavy elements is a real 
process in these stars. The fact that iron-peak elements are 
enhanced below the surface does not mean that the surface 
abundances are necessarily enhanced, as evidenced by the 
observations by O'Toole & Heber (2006). 

The models are robust and predictive: unstable pulsa- 
tions require levitated iron near (but below) the surface at 
the location in the H-R diagram where we find pulsators. 
This permits the blue (and red) edges of the instability re- 
gion to constrain the internal iron distribution in that mod- 
els that are unstable should match (in Toff and log g the ob- 
served instability region. This form of "nonadiabatic aster- 
oseismology" is explored by Charpinet et al. (2008a, 2009), 
who show that models with radiative levitation do indeed 
match the observed instability region for the short-period 
sdB pulsators. 

4 Core properties of sdB stars via seismology 

While much of the activity in studying the seismology of 
sdB stars has been concentrated on the p-mode pulsators 
and diagnostics of global and envelope properties, the long- 
period pulsators hold promise for studying the cores of sdB 
stars. This has been a relatively limited area because of the 
difficulties of observing and resolving dense pulsation spec- 
tra with periods of thousands of seconds (see, for example, 
the heroic effort by Randall et al. 2006), but Kepler obser- 
vations of these stars will provide ample data for g-mode 
probes of their cores. Stay tuned! 

In the interim, those sdB pulsators with periods in the 
300s to 600s range (such as Feige 48, PG 1605H-072, V391 
Peg, and Balloon 090100001) do give us some hope of test- 
ing models of their cores. Measurement of rates of period 
change in these stars can constrain their evolutionary status 
and therefore provide information on their core configura- 
tions. Also, for evolved sdB stars (beyond helium core ex- 
haustion) many of the observed longer-period modes are of 
mixed character - that is, they behave as g-modes in the core 
but p-modes in the envelope. Within their observed period 



M = 0.4780M3, M^„=0.002 Mg 



5.0 




5. 1 


'''expand \ 

L Rj3xl0'y\ 




5.2 






5.3 




j "^expand / 

S32xl0»y/ 


5.4 


'''expand 

- , , , 1 ' 





30000 28000 26000 24000 22000 
T 

'eft 

Fig. 3 Evolutionary track for the same sdB model se- 
quence as in FiglT] The bold portions of the track denote 
times when the radius of the model is increasing. The time 
scale for expansion, i?* /R, is indicated. During core helium 
burning, the radius increases on a time scale that is an order 
of magnitude slower than the expansion rate following the 
establishment of the heUum-burning shell. 

range, there are a mixture of modes that sample very differ- 
ent depths. If we can identify those modes, they can provide 
seismic probes of the core. 

4.1 Rate of period change - secular evolution 

One of the most exciting side-results of sdB seismology 
was the discovery by Silvotti et al. (2007 - hereafter S07) 
of reflex orbital motion induced by a planetary companion. 
While the planet discovery is truly important, the measure- 
ment of the rate of period change (P or dP/ dt) for two in- 
dependent modes in the star is also notable as the first such 
measurement in an sdB star. The values for the two modes 
are shown in Table [T] 

In the short-period sdB stars, since the pulsation fre- 
quencies scale inversely with the square root of the density 
P is determined by the rate of contraction (or expansion): 

In addition to this general scaling, the values of P also de- 
pend on the eigenfunctions. The shape of the eigenfunctions 
determine what parts of the star are most important in deter- 
mining the pulsation frequencies and the rates at which they 
change. For V391 Peg, the values of P differ by 2a. 

For V391 Peg, the time scale for period change {P/ P) is 
roughly H-6.5xl0^y as measured by S07. This corresponds 
to a time scale for radius increase of approximately lO^y, 
which is much faster than the evolutionary time scale for 
sdB stars that are burning helium in the core. Furthermore, it 
is inconsistent with any evolutionary phase where the stellar 
radius is decreasing with time. 

In Fig. |3] we show the same evolutionary track as in 
Fig. [T] with highlighted portions indicating evolution with 
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Table 1 Periods, rates of period change, and model com- 
parisons for V391 Peg. 



Quantity 


Observed 


Shell-burnning 


Core-burning 




(S07) 


114.8 Myr 


88.3 Myr 


T fr 
J- cn 


29300 ± 500 


29219 


28395 


logs 


5.4 ±0.1 


5.358 


5.507 


M, 




O.478OM0 


O.4902M0 


Mh 




O.OO2OM0 


O.OOO2M0 


Pi 


349.54 s 


348.69 s 


349.51 s 




+1.46 ±0.07 


+1.55 


+0.028 


P/P 


7.6xl0®y 


7.2xl0®y 


4.0xl0**y 


{l,n) 




(1,-3) 


(1,1) 


P2 


354.10 s 


352.96 s 


353.14 s 


P(W'^^) 


+2.05 ±0.26 


+2.67 


+0.029 


P/P 


5.5xl0''y 


4.2xl0''y 


3.9xl0*y 


{l,n) 




(0,1) 


(0,0) 



increasing radius, along with the time scale for radius in- 
crease. To match the observed values of P of V391 Peg, 
the models suggest that this star is well beyond the core he- 
lium burning phase, and is evolving towards the white dwarf 
cooling phase with a helium burning shell surrounding the 
helium-exhausted core. 

It is interesting to note that models can be found that 
match the two main periods in V391 Peg, as well as the 
observed T^s and log g, but which are in very different evo- 
lutionary stages. Table [T] shows two representative models 
- one burning helium in a shell (similar to the 114.5 Myr 
model in Figs. [T] and |2|i and one still burning helium in the 
core (similar to the 96 Myr model but with a thinner hydro- 
gen shell). However, to match the observed values of P, the 
shell-burning model is very close while the core-burning 
model is changing too slowly (by a factor of 100). Thus 
measured rates of period change can place important con- 
straints on the evolutionary status of the sdB pulsators. 

While this sketch of the way that P can help constrain 
modes looks promising, it is important to also note that the 
measure rates of period change will always be upper limits 
to the evolutionary rates of the stars. All secular processes 
occurring within the stars operate on the longest available 
time scales, so any process that might influence the phase or 
frequency of a mode will push the resulting determination 
of P to faster changes. Nonlinear mode coupling, stochastic 
phase variations or other "weather-like" phenomena could 
make it appear that the star is evolving more quickly. There- 
fore, the conclusion that V319 Peg is highly evolved re- 
quires more precise modeling and (hopefully) observations 
of P for additional modes. 

4.2 Mixed-character modes 

While the short-period sdB variables are generally consid- 
ered to be p-mode pulsators, with the oscillations being in- 
sensitive to core properties, the lower-gravity members of 
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Fig. 4 Properties of mixed-character modes in an evolved 
sdB model that has exhausted helium in its core. Top: 
the propagation diagram, with the Brunt-Vaisala frequency 
(solid line) and Lamb frequency (dashed line) as a function 
of fractional radius. The horizontal line Ues at the frequency 
corresponding to the ^=1, 7i=-3 g-mode; diamonds denote 
nodes in the radial eigenfunction that are g-mode-like, and 
the circle shows the envelope node in the p-mode propaga- 
tion region. Bottom: the integrated adiabatic period weight 
functions. For the /=0, n=\ radial mode, the node lies at 
r jR = 0.77. The period of pulsation is 348.7s for the ^=1 
mode, and 353.0s for the Z=0 mode. In both panels, the dot- 
ted line shows the hehum composition (with vertical scale 
from to 1). 

this class generally show longer-period pulsations (i.e. in 
excess of 300 s). For these stars, the models indicate that 
some of the modes are mixed-character, as has been shown 
by Kawaler (1999) and, more extensively, by Charpinet et 
al. (2000) and Hu et al. (2009). Some modes have oscilla- 
tion eigenfunctions have nodes close to the center and be- 
have like 5— modes there, but also have a node or nodes 
near the surface, where the eigenfunction is more p-mode- 
like. That means that some the oscillations in these stars can 
allow us to probe the core. 

Fig.lHshows some internal properties of a model of such 
a star - the model that provides a match to the period and P 
of V391 Peg listed in Table [T] In the top panel, the Brunt- 
Vaisala frequency shows two bumps at the position of the 
H/He transition and at the edge of the He-exhausted core. 
For the Z=l, n=-3 mode, the inner composition transition 
acts as a reflecting boundary, partially trapping a mode in 
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Fig. 5 Avoided crossings in an evolving, helium-shell 
burning sdB model. Period is shown vs. log g, and the model 
evolves from right to left. Horizontal Unes indicate periods 
of the two main modes in V391 Peg. Each vertical string 
of dots represents periods in a model at a given log g; lines 
connect modes with same values of n and I. The complex 
behavior is characteristic of avoided crossings. The vertical 
(blue) dashed line is the value of log g that has the closes 
period match between the observations and model; this is 
the model detailed in the third column of Table [T] 

the interior, and enhancing the period weight function sig- 
nificantly in that region. Even though that mode has a fre- 
quency that is very similar to the frequency of the first over- 
tone radial mode, it is in essence a g-mode. 

The existence of mixed-character modes in sdB stars 
leads to avoided crossings and other phenomena associated 
with mode trapping. Fig.|5]shows this complex behavior in 
the shell-burning O.478M0 model. This figure clearly shows 
avoided crossings. While this complicates the mode identi- 
fication and analysis, the value of mixed character modes 
in probing the deep internal structure of sdB stars has enor- 
mous potential. In addition, it provides an opportunity to 
examine the rotational state of the core to constrain mod- 
els of angular momentum transport on the HB and during 
earlier stages (Kawaler & Hostler 2005). 

5 Conclusions 

Through asteroseismology, we can begin to answer basic 
questions about stars like our own Sun that have perplexed 
stellar astrophysics for four decades. Not only is the origin 
of the hot subdwarf B stars one of the remaining myster- 
ies of low-mass stellar evolution, but the stars themselves 
are unique laboratories for studying the physics of diffusion, 
radiative levitation, hot stellar winds, turbulence, core con- 
vection, and semiconvection. Studies of the cores of these 
stars can provide a window to these important physical in- 
gredients of a wide variety of stars ranging across the entire 
H-R diagram. 

Many of these stars are in binary systems, and in fact 
binary evolution may be required to understand their origin. 



Since many do reside in close binaries where tidal effects 
are substantial, we can use the pulsations to probe processes 
such as interactions between pulsations and tides (Reed et 
al. 2006), tidal circularization, and spin-orbit synchroniza- 
tion (see for example Charpinet et al. 2008b, van Grootel et 
al. 2008). 

As exciting as the developments in sdB asteroseismol- 
ogy have been, we anticipate that space-based observations 
with Kepler will provide revolutionary new data on pulsat- 
ing sdB stars - in particular, the long-period g-mode pul- 
sators. With these observations in hand, we can look for- 
ward to peering at never-before-seen details of the deep in- 
terior of stars that represent the future fate of our own Sun. 
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